Gravitational Time Delay Effects on CMB Anisotropics 
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We study the effect of gravitational time delay on the power spectra and bispectra of the cosmic 
microwave background (CMB) temperature and polarization anisotropies. The time delay effect 
modulates the spatial surface at recombination on which temperature anisotropies are observed, 
typically by ~ 1 Mpc. While this is a relatively large shift, its observable effects in the temperature 
and polarization fields are suppressed by geometric considerations. The leading order effect is from 
its correlation with the closely related gravitational lensing effect. The change to the temperature- 
polarization cross power spectrum is of order 0.1% and is hence comparable to the cosmic variance 
for the power in the multipoles around £ ~ 1000. While unlikely to be extracted from the data in 
its own right, its omission in modeling would produce a systematic error comparable to this limiting 
statistical error and, in principle, is relevant for future high precision experiments. Contributions to 
the bispectra result mainly from correlations with the Sachs- Wolfe effect and may safely be neglected 
in a low density universe. 
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I. INTRODUCTION 

In order that the full potential of anisotropies in the 
cosmic microwave background (CMB) temperature and 
polarization fields be realized, effects that have been pre- 
viously dismissed as negligible in their own right must 
now be reconsidered as potential sources of systematic 
error. Projections as to the ability of CMB experiments 
to measure fundamental cosmological quantities jl] pre- 
cisely and reveal information about the structure forma- 
tion process from secondary effects Q rely on the fact 
that statistical errors from the sampling of a finite sky 
rapidly decrease toward smaller angular scales. Statis- 
tical errors in the power spectra decline from ~ 1% at 
degree scales to ^ 0.1% at the several arcminute scale. 
To achieve this precision in practice, all physical, astro- 
physical and instrumental effects at this level must be 
included in the analysis to avoid generating systematic 
errors that are comparable to the statistical errors. 

A host of physical effects contribute to the anisotropies 
at second order in perturbation theory 1^-^. Since pri- 
mary anisotropies are formed at recombination when the 
cosmological density perturbations are at the 10~^ level, 
most second order effects are entirely negligible. There 
are two general ways in which higher order effects can be 
important. Firstly, the primordial perturbations respon- 
sible for the primary anisotropies grow into non-linear 
structures today by gravitational instability. Effects that 
take advantage of this fact mainly involve scattering of 
CMB photons at low redshifts in large-scale structure 
and non-linear objects [0-0. Secondly, since recombi- 
nation, CMB photons propagate across essentially the 
whole horizon volume. Intrinsically small effects can ac- 
cumulate along the path. Indeed it is well known that the 
gravitational lensing of CMB photons has a substantial 
effect on the power spectrum of the anisotropies [|ll|,|l^ . 

In addition to the lensing effect, gravitational poten- 
tials of large-scale structure contributes a time-delay ||l^ 



that accumulates along the path - an effect familiar from 
studies of the light-curves of lensed quasars (e.g. iQ). 
In the case of the CMB, the time-delay warps the spa- 
tial surface at recombination from which the primary 
anisotropies arise Because the lensing depends on 
the angular gradient of the projected potentials whereas 
the delay depends on the projected potential itself, the 
fractional delay is generically smaller than lensing and 
has not been explicitly calculated in the literature. We 
shall see however that because of the angular smoothness 
(or coherence) of the lensing, the reduction in amplitude 
is not in and of itself large. Furthermore, gravitational 
time delays are strongly correlated with lensing, lead- 
ing to additional effects in the power spectrum. Indeed, 
the typical perturbation in comoving units is on the or- 
der of 1 Mpc. The effect of gravitational time-delay on 
the spectra of temperature and polarization anisotropies 
therefore merits further study. 

In §^we present the formalism required to understand 
these gravitational effects on the temperature and polar- 
ization fields. We proceed in §111 and |^ to evaluate the 
delay and lensing-delay correlation effects on the power 
spectra of tem pe rature and polarization anisotropies re- 
spectively. In §M, we consider their effects on the bispec- 



tra (three point correlations). We conclude in §|VI| with 
a discussion of our results. 

To illustrate our calculations, we assume a cold dark 
matter model (CDM) with a cosmological constant (A) 
with parameters flc = 0.30 for the CDM density, fib = 
0.05 for the baryon density, Oa = 0.65 for the vacuum 
density, h — 0.65 for the dimensionless Hubble constant 
and a scale invariant spectrum of primordial fluctuations, 
normalized to COBE llTl. 



1 



^ecombinafTo^ 



source 




FIG. 1. Gravitational Lensing vs. Time Delay. Lans- 
ing introduces an angular perturbation in the mapping of a 
plane-wave source field at recombination onto anisotropies to- 
day. Time delay introduces a radial modulation. When the 
wavevector is perpendicular to the line-of-sight, features in the 
angular spectrum - such as the acoustic peaks - are created, 
geometrically distinguishing the otherwise similar lensing and 
delay effects. 



II. FORMALISM 
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FIG. 2. Power spectra for the lensing deflection angles 
{(f>4>), time-delay (dd) and deflection-delay cross correlation 
{(pd). The underlying lensing potential spectrum C^** and 
cross spectrum Cf"^ are weighted by £{£ + 1) and [£(£ + f )]^/^ 
respectively to reflect the angular gradients in the deflection 
angles. 



As the CMB photons propagate to the observer from 
the recombination epoch {z ~ 10"^) through large-scale 
structure in the universe, they suffer the effects of grav- 
itational lensing and time delay. These effects are both 
formally second order in perturbation theory because 
they would leave a homogeneous and isotropic CMB un- 
perturbed (c.f. U&). There are a host of other second 
order effects |^-||^ To the extent that they are uncorre- 
lated with each other, they may be viewed as independent 
effects. The lensing and time delay effects are strongly 
correlated because both arise from the gravitational po- 
tentials of large scale structure and must be considered 
together. We therefore begin with a review of lensing 
effects. 



A. Lensing 

Lensing involves a deflection that remaps the tem- 
perature and polarization fields according to the angu- 
lar gradient of the lensing-weighted projected potential 
n ^ n V0(n) (see Fig. 

0(n) = -2y" d??50('7)*(?-n,r?), (1) 

where 

I rv r' -r 

/ dri're-^^. (2) 
'"Jo ^ 

Here overdots represent derivatives with respect to con- 
formal time — J dt/a and <i>(x, 77) is the gravitational 
potential. In an open universe, the conformal distance 
traveled by a photon rij]) = rjo — rj should be replaced 



by angular diameter distances. We take comoving units 
where c — 1 throughout. The projected potential itself is 
a field on the sky and may be decomposed into multipole 
moments 

0(n) = ^0,„yr(n), (3) 

£m 

and described by its power spectrum 

{(f'im't'i'rn') = SiJ'Sm,m'Cf'^ . (4) 

This power spectrum is shown in Fig. ^ 

A useful measure of the amplitude of the lensing effects 
is the rms deflection angle 6^^^ as defined by 

e=i 

The factor £{£+!) reflects the angular gradient in the def- 
inition of the deflection. Note that the £ = monopole 
does not contribute since its angular gradient vanishes. 
In the fiducial ACDM model e,n,s = 7.5 x lO"'^ or 2.6'. 
Note that this is in sharp contrast with the angular co- 
herence of the deflection angle. The variance 0^^^ reaches 
half its total value by £1/2 = 30 in the fiducial model or 
an angle of 6*1/2 = '^7^/^1/2 ~ 0.02 (1°). The smaller 
scale potential fluctuations tend to produce deflections 
that cancel out along the line of sight. 

B. Time Delay 

Now consider the time-delay effect. Photons follow null 
geodesies in the perturbed metric so that in the fixed time 



2 



interval since last scattering, the distance traveled by the 
photons is perturbed as r ^ r[l + d(n)] (see Fig. |^) where 



2 [ 

(i(n) = / d'qe '^$(rn, 

J 



(6) 



This is referred to in the literature as the potential or 
Shapiro time delay. The geometric time delay is of order 
the deflection angle squared and is hence substantially 
smaller for time-delays across angular scales much larger 
than 6'i.ms- 

The delay field may also be expanded in spherical har- 
monics 



and characterized by a power spectrum 



dd 



(7) 



(8) 



Although both the lensing and time-delay effects are 
based on the gravitational potential projected along the 
line of sight, there is an important difference between 
the two. Lensing depends on the angular gradient of 
the potential and hence its observable consequences are 
weighted by 1). This has the effect of increasing 

the magnitude of the effects and weighting it to higher 
multipoles. 

We can see these effects in the rms delay llSl 



"rms / J 



21+1 
in 



-c; 



dd 



(9) 



For the power spectrum of the fiducial ACDM model (see 
Fig. I) 4nis = 5.4x10"^ or - 0.5/1"^ Mpc. Although this 
is a relatively large shift, we shall see that its observable 
consequences are reduced due to the large coherence scale 



of the delay. 
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Given these considerations as to the magnitude and co- 
herence scale of the effects, delay effects can be enhanced 
through their cross-correlation with lensing 



(10) 



The two fields tend to be well correlated as can be seen 
by the rms 



(11) 



where Crms = 1-2 x 10 and £1/2 — 5 in the fiducial 
ACDM model (see Fig. |). To evaluate the observable 
consequences, we now turn to the angular and spatial 
structure of the temperature and polarization fields. 
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FIG. 3. Projection functions for the lensing and delay ef- 
fects. Lensing involves je{kr) whose strong peak in £ can 
retain source features; the j'likr) of the time-delay cannot. 
The product of these functions reflects the cross correlation 
and is suppressed due to the phase difi'erence. 



C. Temperature Field 

The CMB temperature field on the sky may be written 
implicitly as the projection of sources S which contribute 
in the optically thin regime and are so weighted by e""^ 
where r is the optical depth (see Fig. |^). In general, these 
sources have intrinsic angular structure of their own and 
are characterized by the spherical harmonic moments of 
their Fourier amplitude S'™'(fc). Explicit forms for the 
sources are given in |17|. 



The contribution from a given wavenumber k to the 
temperature field on the sky today may be formally ex- 
pressed as 



(■no 

e(n; k) = / di^e-^ ^ 5™' (77; fc)G™' (rn; k) , (12 



where r \ 



Gr(x;k) = (-*) 



2£+l 



y/"(n)cxp(ik-x). (13) 



In an open universe, the plane waves must be replaced 
with eigenfunctions of the Laplacian in a curved space. 
Note that we will often omit the /c-index where no con- 
fusion will arise. 

The angular structure of these relations can be simpli- 
fied by considering a specific frame where z || k and 

exp(ik • x) = ^(-i)V47r(2^ + l)j£(fcr)Yg"(n) . (14) 

Provided that the angular basis does not change in tran- 
sit, one can then sum up the orbital angular momentum 
from the plane wave with the intrinsic angular momen- 
tum of the source. This assumption is the angular equiv- 
alent of the "Born approximation" where the lensing is 
evaluated on unperturbed trajectories; we shall see in 
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the next section that it is a good approximation for the 
lensing and time-delay effects on sources with low order 
intrinsic angular structure. 

In this special basis, the product of the intrinsic (i^™') 
and plane wave (Y^) angular momentum may be reex- 
pressed through the addition of angular momentum. The 
temperature field then becomes flTIl 



e(n;k) = ^ef„(fc)yr(n) 



(15) 



Here, the operator 

^mi [jl] = 



(16) 



and j^'™' are linear combinations of with weights given 
by the Clebsch-Gordon coefficients of the coupling. The 
fundamental functions are = ji (for isotropic pertur- 
bations, e.g. gravitational potentials) and 



,■22 



l3{e + 2y. 1 

8 (^- 2)! 2^-1-1 



je-2 + je je + je+2 
(2£-l) (2£ + 3) 



, (17) 



(for transverse quadrupole sources, e.g. gravitational 
waves). Note j^^ oc ji{x)/x'^] we have written it out 
here to emphasize linearity. Others are given in [0 , but 
can be written in terms of these fundamental functions 
through integration by parts [22| . 

Since the basis for the expansion is linked to the di- 
rection of k, integrating over modes to obtain the final 
power spectrum in principle requires a series of rotations 
into a fixed basis. In practice, the statistical homogene- 
ity of the source and isotropy of the angular distribution 
requires that the fc-modes add in quadrature and the an- 
gular power spectrum is independent of m so that we 
may replace individual multipoles with an average over 

TO 



dk 



1 



k 27r2 2^ -f 1 ^ 



(e;„(fc)e,„(fc)) . (is) 



For the primary anisotropics, the power spectrum is sim- 
ply Cf^ = T^oo^ ^YiCTC 



rpab 

1 (/ — 



dk fc3 (fc,yo)-+^ 

I 2^ ^nii I 



k 27r2 2e- 



Im, [j. 



(6)1 



(19) 



Here we have used the short hand convention that j^"^ = 
jf, j^^-* = and j'f^ = j" with primes denoting deriva- 
tives with respect to the argument. 

For the higher order lensing and delay effects, we ex- 
pand equation ( p^ ) for the temperature field to second 
order in the relevant quantities 
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FIG. 4. Power spectra for the angular {e{£+l)Ti°), radial 



(r/^) and angular-radial cross gradients {^/I(I+l)T^'^) of 
the primary anisotropics. The cross gradient spectrum has 
been weighted to reflect the I contributions from the angular 
gradient. 



e(n) = e"(n) + e'^(n) + e'^(n) 

+ ^e*' (n) + e*'^(n) + ie'^' (n) , (20) 

where 6*^(11) is the zeroth order contribution from the 
primary anisotropics. 



(21) 



e^(n) 


errii 


[je] V,0(n)V^rr'(n), 


e''(n) 


= E^™ 

errii 


[f,]{k7^o)din)Yrin), ( 


e^'(n) 


= E^™ 

errii 


[je] V,0(n) Vj0(n) W^Yp{h), 


e^''(n) 


= E^™ 

errii 


[f,] (A:%)d(n)V.0(n)V^yr(A), 


e'*'(n) 




[f;] {krjordin)d{h)Ypin). 



errii 

We shall see that the evaluation of these effects reduces 
to the computation of the higher order derivative power 
spectra in equation (|l9[). These power spectra are easily 
evaluated with minimal modifications to the publically 
available CMBFAST code: we simply replace ji with j'^ 
and j" and leave the evolution and integration over the 
sources unchanged. To the extent that the underlying 
sources themselves are smooth, integration by parts on 
equation (|l^ shows that t/"^^^^''^^' « e.g. T/^ « 

—T^"^ = —T^^. This also implies that terms such as 
= are suppressed; mathematically this is due 
to the lack of correlation between ji and j'p (see Fig. H). 
These spectra are compared in Fig. |4| 

There is an additional effect from lensing due to the 
fact that sources are in general anisotropic and gravita- 
tional lensing changes the angle at which they are viewed. 
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This formally violates the assumption that the orbital 
and intrinsic angular momentum of sources can simply 
be added without additional remapping. The anisotropy 
in the actual sources, however, are confined to ^ < 2. To 
obtain an order unity effect, the deflections must change 
the viewing angle by order unity in radians. To obtain an 
order unity effect for lensing, the deflections need change 
the viewing direction by only a wavelength of the angular 
perturbation £~^. Thus at sufficiently high observed £, 
deflections always win. 



D. Polarization Field 

The complex Stokes parameter P± = Q zL ill of the 
polarization on the sky can be expressed implicitly as 
the projection of the quadrupole moments of the photon 
temperature and polarization distributions at last scat- 
tering 

P±(n;k)=/ dr;e--^Q'"'(,7;fc)±2Gr(rn;k), (22) 
Jo 



where Q"^- = -f(V6e2m - 6E2m)/lO (see 0, for an 
explicit definition) and 



±2Gr(x;k) = (-*)' 



An 



2i+l 



±2r2'"(n)exp(zk-x). (23) 



Here, ±2^7" spin-2 spherical harmonics |]l9[. A 

recoupling of the spin spherical harmonics as in equa- 
tion (nl) yields 



P±(n;k) = J2iEf-n{k) ± ^Be„^{k)] ±2Yr{n) 

irn 

= ^P„J±a,]±2yr(n), (24) 

enii 

where the operator 

rvo 

x^Q™'(,?;fe)±ar(M- (25) 

Here, ±a™' — e™' ± where the latter are combi- 
nations of ji given in p7| that define the projection of 
the source onto the E and B polarization modes. For 
quadrupole sources related to density (mi = 0) fluctua- 
tions, e° = jj^ (see eqn. |l^) and /3° = so that there 
are no contributions to P-parity polarization. 

The power spectra for the polarization and temp- 
erature-polarization cross correlation are defined as 

=Jt&2iTi^ {FLik)FUk)) . (26) 



where F and F' take on the values Q,E and B. For the 
primary polarization, we define Cf^ = P"°, Cf^ = P^°, 
Cf ^ = X^oo, where 



Ef = 



dk P (fcryo)'' 



k 27r2 2£ 



ab 



B 



dk P (fc^„)a+b 



dk P (fcryo)'^ 



k 27r2 2£- 



Other combinations vanish due to parity considerations. 
As in equation ( |l9|) the indices a and b refer to deriva- 
tives of the underlying Bcsscl functions. The higher or- 
der gravitational effects on the polarization field may be 
evaluated through a second order expansion of the polar- 
ization field 

P±(n)=pO(n)+P|(n) + P|(n) 

+ ipf (n) + pf(n) + ipf (n), (28) 

where P±(n) is the zeroth order contribution from the 
primary anisotropics, 



P|(n) 






V,^(n)V^±2i^r(ft), 


(29) 


Plih) 


= E^- 

Inii 




(fcr;o)d(n)±2r/''(n). 




Pf in) 


= E^- 

errii 


[±ae_ 


V,0(n) V,0(n)V''VJ±2rr 


'(ft), 


Pf{h) 


~ E ^'"i 
inii 




(fcr;o)d(n)V,(/)(n) V'±2yr 


(n). 


pf (n) 


= E ^'"i 




(A:%)2d(n)d(n)±2l^,"'(n) 





The evaluation of the effects will involve the higher or- 
der derivative power spectra of equation (|2^). Again we 
modify CMBFAST to calculate these spectra. Just as for 
the temperature power spectra jff-^^^''^^^^'^ ~ —E^^ and 
similarly for B and X. This also implies that terms such 
as = are suppressed. 



III. TEMPERATURE POWER SPECTRUM 

The perturbations to the power spectrum due to the 
lensing and time-delay effects follow by considering the 
second order terms in the two-point correlation of the 
temperature field in equation (pO|). We can express the 
contributions schematically as 

cf Q - T,"" + (t/^ + pf ) + (T/'^ + rf ) + t/'^ . (30) 

We will now define and consider each term in turn. 
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FIG. 5. Delay, lensing and delay-lensing (cross) perturba- 
tions to the CMB temperature power spectrum for the fidu- 
cial ACDM model. The cross spectrum dominates the delay 
spectrum but both produce negligible changes to the primary 
anisotropics unlike lensing. 



A. Lensing Spectra 

The pure lensing contributions can be defined as 
dk 



rn 

k 27r2 2^ -f 1 ^ 2 V\ 



t: 



(31) 



Following Ref. [gO| , we expand the perturbations to the 
temperature field fi'om equation ( pT[ ) and the projected 
potential (j) in spherical harmonics to obtain 



rp4"t> _ \ " r"^0rT^OOTT/2 r 2 



(32) 



where 



(2^1 -fl) (2^2 + 1) 



47r 



I h h 




(33) 



The Weeii2 term comes from the integral over the product 
of three spherical harmonics; the Luit2 term comes from 
the conversion of angular gradients into angular Lapla- 
cians through integration by parts. On scales where there 
is power in the primary power spectrum, the two terms in 
equation ( |3^ nearly cancel. The reason is that the large- 
angle modulation produced by lensing simply smoothes 
the features in the power spectrum across A£ ^j^yg ~ 30 
[ pl| . The result of combining the two terms in the fiducial 
ACDM model is displayed in Fig. |. 



B. Delay Spectra 

The time-delay modifications to the power spectra can 
be defined as 

dk k^ 



m 



dk k^ 
T2 



(34) 



and their evaluation follows closely that of lensing except 
for the simplification that no derivatives of spherical har- 
monics are involved. The result is 



eii2 



rpd^ t2 rp{)2 

where we have used the identity 

2h + 1 



ETAmi*-t^mi 



to evaluate 



E 



dhYr*Yr'*Yr'Y„ 



2^1 + 1 

( 

4tt 



(35) 



(36) 



(37) 



The main difference between the lensing and time delay 
contributions is that the power spectra of the angular 
gradients of the projected potential £{i + l)Cf'^ and pri- 
mary temperature anisotropies £(£ + 1)T™ are replaced 
by the power spectra of the delay C^*^ and that of the 
radial derivatives of the temperature field T^^ « . 

As in the case of lensing, the time-delay effect repre- 
sents a smoothing of the gradient power spectrum T/^ 
across a width of A£ ^ £ij2 ~ 2. The effect is strongly 
suppressed by the small width of the smoothing and the 
nearly featureless underlying spectrum Tj^ (see Fig. Q). 
Contrast this with lensing which has a smoothing width 
of ~ 30 and an angular gradient power spectrum that 
shows peaks from the acoustic oscillations. Acoustic fea- 
tures in the power spectrum arise from features in the 
source power spectrum when the source wavevector is 
oriented perpendicular to the line of sight. The source 
then projects with a one-to-one correspondence between 
angular and physical scale {£ = kr see Fig. |l|). Other 
alignments contribute to a broad tail to lower multipole 
moments (see Fig. In the perpendicular orientation, 
however, a perturbation to the radial distance simply 
moves the last scattering surface along the crests and 
troughs of the source leaving no net effect. Mathemati- 
cally, this effect can be seen in the fact that ji{kr) as a 
function of £ possesses a strong peak at £ = kr, whereas 
j'i{kr) does not (see Fig. ||). 

The net effect is therefore much smaller than the naive 
scaling of the lensing rms ^ and delay rms (||) would 
imply. It is shown in Fig. ^ for the fiducial model. 
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C. Cross Spectra 

The cross correlation between the lensing and delay 
cause modifications defined by 

rp<i>d ^ f dk 1 ^ /- 



_ / ClP* p) 

~2 K^lm^irn 



CC, 



(38) 



The e'^'^ terms are identically zero since 



= 0, 



(39) 



by virtue of the addition theorem of spherical harmonics. 
The first term reduces to 



'ct''T?}Wl,,L 



(40) 



Unlike the pure lensing and time delay contributions, 
there is no second canceling term. However as discussed 
in § p C| , T^"^ is intrinsically small reflecting the lack of 
correlation between and j[ (see Fig. Nonethe- 
less, its net contribution to temperature anisotropy power 
spectrum is still larger than the pure delay contribution, 
as a result of the larger amplitude and smaller coherence 
of the lensing-delay power spectrum Cf'^ (see Fig. ||). 



IV. POLARIZATION POWER SPECTRA 

Evaluation of the lensing and delay effects on the polar- 
ization and temperature-polarization cross power spectra 
follows precisely the same steps as that of the tempera- 
ture power spectrum considered above. We can express 
the contributions schematically as 



r<BB nOO 

riBE _ vOO 
— 



Et)^{Ef + Ef)^ 
Bf) + {Bf + Bf)^ 

-xf) + {xr^xf) 



vxf 



(41) 



We will now consider each term in turn. 



A. Lensing Spectra 

Following the lensing contributions to the polar- 
ization power spectra can be evaluated as 



1 




(42) 
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FIG. 6. Delay, lensing and delay-lensing (cross) contribu- 
tions to the CMB polarization power spectra for the fiducial 
ACDM model. The delay-lensing cross contributions to the 
temperature-polarization power spectrum reaches 10^^ of the 
primary and/or lensing power. Delay effects only weakly gen- 
erate power in the BB spectrum. 



f„,(^2^^-4)S,, (43) 



<;i<;2 



7 



X, 



1 



where 



± 
Here 



'-ET ± 



(2£i + l)(2^2 + l) 



47r 



4 ^2 

-2 



(44) 
(45) 
(46) 



comes from the integral over the product of a spherical 
harmonic with 2 spin-2 spherical harmonics. As in the 
case of the temperature power spectra, the main effect 
on the QE and EE power spectra is a smoothing by 



M 



If there were an intrinsic BB power spec- 



trum, the smoothing would also apply. However, since 
scalar perturbations do not generate BB modes in the 
primary polarization {Bf' = 0), the generation of B- 
polarization from the lensing modulation of the primary 
-E-polarization dominates |12| . The amount of genera- 
tion is again related to the coherence scale lij2 ~ (-i of 
the effect as reflected in the difference between even and 
odd L terms in Vu-^i^. Considering the triplet as forming 
a triangle, even terms are associated with the cosine of 
twice the opening angle between £ and £2; odd terms are 
associated with the sine of that angle ( |2^, Eqns. (B8) 
and (BIO)). These polarization contributions are shown 
in Fig. 6 for the fiducial model. 



B. Delay Spectra 

The derivation of the delay power spectra follows the 
same steps as those involved in the lensing derivation 
yielding 



Ef 



drmsE' 



02 



B? 



02 



vdd 

JLf — 



XI 



02 



(47) 



(48) 



(49) 



The results of summing these nearly canceling pairs for 
the fiducial model are shown in Fig. 6. A i?S-spectrum 
is generated out of the primary E'.E-spectrum but at an 
efficiency that is substantially below that of gravitational 
lensing. The underlying reason again is that the coher- 
ence of the effect £1/2 ^ 2. 
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FIG. 7. Delay-temperature cross correlation. The 
cross-correlation between the temperature and delay fields, 
as relevant for the bispectra, are dominated by contributions 
from the Sachs- Wolfe (SW) efi'ect with partially canceling con- 
tributions from the integrated Sachs- Wolfe effect (ISW). 



C. Cross Spectra 

The cross spectra contributions are defined as 
Ef" = 2Y,Cf^^Y.f^{£ -f 4 + mIi,iMt,i2 , 

Bf = 2Y,Ct^^Y.l{t + £,+ ^2)V,1,,L«,,, , 
I1I2 

Xf = ^t'iX^^l + xr^ )Wu,,,Vu,i2Lu,i2 ■ (50) 



Figure 6 shows that the cross spectra dominate over the 
pure delay spectra for the QE and EE power spectra. 
In particular, since the lensing effects themselves ap- 
proach order unity at ^ '--^ 1000, the lensing-delay cross 
effects reach ~ 10^'^ of the primary Cf^ power spec- 
trum. While still small, the contribution is of order the 
cosmic variance out to comparable multipoles. 



V. DELAY BISPECTRUM 

Second order effects generally produce non-Gaus- 
sianity in the CMB temperature and polarization fields. 
Effects that provide a negligible change to the power 
spectrum can in principle produce observable effects due 
to the expected Gaussianity of the primary anisotropics. 
Here we consider the three-point correlations induced by 
time-delay in angular harmonic space, i.e. the bispec- 
trum. 

The angle averaged bispectrum is defined as 
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FIG. 8. Signal-to-noise for detection of bispectra involving 
gravitational time-delay in an ideal cosmic-variance limited 
experiment. 
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(51) 



where the F's can take on the values O, E, B for the 
temperature and polarization components respectively. 
Since the derivation follows closely that of the lensing 
bispectra terms, we refer the reader to ||2^ for detailed 
derivations. 



A. Temperature 

The temperature bispectrum produced by time de- 
lays follows immediately from the decomposition in equa- 
tion (|2ll) 



B\ 



eee 



v/2£i + W^,fe^3^^W' + 5perm. , (52) 



where the permutations are with respect to the ^-indices 
and 



Bd 



(53) 



Unlike lensing however, Cf"^ is not solely the result of 
correlations with secondary anisotropics but includes rel- 
atively large contributions from the primary anisotropics 
themselves. The delay field arises from potential fluctu- 
ations of sufficiently large scale that it is correlated with 
the Sachs- Wolfe effect in the primary anisotropics. 
To see this, let us approximate the primary anisotropics 
at large angles as the sum of Sachs- Wolfe (SW) and in- 
tegrated Sachs- Wolfe (ISW) contributions 



e(n) 



1 



$(r,n, ?7,) - / dr]' 2$(r'n, ?/) , (54) 



where the asterisk denotes evaluation at recombination. 
This temperature field may be directly correlated with 
the delay field in equation (^) . The resulting power spec- 
trum is shown in Fig. |. Note that the SW and ISW 
contributions cancel at the lowest multipoles. Unlike the 
lensing-ISW correlation 1^, the delay effects on the bis- 
pectrum do not vanish and indeed increase as J^a 0. 

To determine whether the contributions are detectable, 
we evaluate the signal-to-noise ratio for an ideal cosmic 
variance limited experiment 



(55) 



We show the cumulative signal-to-noise out to a given 
maximum £ in Fig. The time-delay temperature bis- 
pectrum is not detectable even for an ideal experiment. 
This should be compared with the lensing-temperature 
bispectrum, where the cross-correlation between lensing 
and effects such as SZ can be detected with upcoming 
experiments, such as Planck. 



B. Polarization 

Bispectrum terms involving the polarization follow 
similarly. The non-vanishing contributions are 



B 



see 



+ ^2£2 + lWe,e,e,X^^C^/ + {£2 <^ £3) , 

'£2), (57) 



^Fefet = ^2^1 + lVi^i^i^Ej°Cfj^ + {£1 



for £1 + £2 + = even and 



see 



= iVW+lVi.e^is Gt^^t + (^2 ^ 4) , (58) 



BEQ 
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Qd 

a ' 



(59) 



for £1 + £2 + £d. ~ odd. The signal-to-noise ratio in an 
ideal, cosmic variance limited experiment is 
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(60) 
(61) 
(62) 
(63) 



In Fig. H, we show the cumulative signal-to-noise ratio as 
a function of the maximum £. The signal-to-noise ratio 
for the BOO term is larger than the others since we have 
assumed that Cf^ vanishes for the primary anisotropics 
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and is only generated by the lensing and delay effects. 
Nonetheless, the signal-to-noise is substantially less than 
unity and so the time-delay effects are unlikely to be 
detectable in the bispectrum or affect the extraction of 
other effects from the bispectrum. 

VI. DISCUSSION 

Gravitational time delays introduce a radial pertur- 
bation in the mapping of the CMB temperature field 
at recombination onto temperature and polarization 
anisotropics today. The effect is closely related to grav- 
itational lensing which introduces angular perturbations 
in the same mapping. Despite the fact that radial pertur- 
bations are only one order of magnitude smaller than an- 
gular perturbations, and moreover highly correlated with 
the angular perturbations, their effect on the power spec- 
tra and bispectra of the temperature and polarization 
fields is approximately 3 orders of magnitudes smaller. 
The underlying reason is that on the angular scales of 
the acoustic peaks neither effect actually generates new 
anisotropics; both induce a large scale modulation of the 
primary field. Angular modulations produce a substan- 
tial effect due to angular structure in the acoustic peaks 
by smoothing the power spectrum on the coherence scale 
A£ ~ 30. Radial modulations produce much smaller ef- 
fects due to the lack of radial structure in the perturba- 
tions that form the acoustic peaks. They also suffer from 
the fact that the angular coherence or smoothing scale of 
the radial modulation is typically one order of magnitude 
larger A£ ~ 2. 

As a result, at ^ ~ 1000, the delay and lensing-delay 
correlation effects are < 10^'* of the temperature, E- 
polarization and lensing-induced B-polarization power 
spectrum generated by lensing. For the temperature- 
polarization cross correlation, the lensing-delay corre- 
lation effect approaches 10"'^ and so is comparable to 
the cosmic variance on these scales. This enhancement 
reflects the same efficiency with which lensing modula- 
tions affect the temperature-polarization correlation. It 
is therefore relevant in principle for the Planck satellite 
or any future experiment that expects to be cosmic vari- 
ance limited at £ > 1000. In practice, there may well be 
other more limiting sources of systematic errors such as 
galactic and extragalactic foregrounds. 

For the bispectra, the time delay couples mainly with 
the Sachs- Wolfe effect in the primary anisotropics. For 
a cosmic variance limited experiment, the signal-to- 
noise is highest in the i?-temperature-temperature bis- 
pectrum since the i?-polarization vanishes for primary 
anisotropics from scalar perturbations. In a realistic ex- 
periment, the low level of the signal will make the cos- 
mic variance limit difficult to achieve. In any case, the 
signal-to-noise ratios in the bispectra never exceed the 
10~^ level for £ 1000 and hence are unlikely to inter- 
fere with the extraction of signals in the bispectrum from 



secondary anisotropics by the next generation of satellites 

The potential of cosmic microwave background 
anisotropics for studying cosmology is considered vast 
primarily because of the physical processes underlying 
their formation are thought to be understood to extraor- 
dinary precision relative to other astrophysical systems. 
Though unlikely to affect the next generation of experi- 
ments, small effects such as the gravitational time delay 
considered here must be calculated and included to en- 
sure that this potential is realized. 
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